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1 Introduction 

The complex interactions existing between turbulence and mixing in a bluff- 
body stabilised flame configuration is investigated by means of a wavelet cross- 
correlation analysis on Large Eddy Simulations. The combined approach allows 
to better point out typical features of unsteady turbulent flows with mixing 
through the characterisation of the processes involved both in time and scales. 
The wavelet cross-correlation analysis of the time signals of velocity and mixture 
fraction fluctuations can be an an effective tool to study the processes involved 
in turbulent mixing flows which are of great interest in combustion problems. 

2 Generalities on wavelet cross-correlation 

The continuous wavelet transform of a function f(t) is defined as the convolution 
between / and a dilated function ijj called wavelet mother: 



where a is the dilation parameter, which plays the same role as the frequency 
in Fourier analysis, and t indicates the translation parameter corresponding to 
the position of the wavelet in the physical space. In the present study we use 
the complex Morlet wavelet (ip(t) = e M °*e _t / 2 ) as wavelet mother. 
Let Wf(a, t) and W g (a, t) be the continuous wavelet transforms of f(t) and g(t). 
We define the wavelet cross-scalogram as 




(1) 



W fg (a,r)=W}(a,r)W g (a,r) 7 



(2) 
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where the symbol * indicates the complex conjugate. When the wavelet mother 
is complex, the wavelet cross-scalogram Wf g (a,r) is also complex and can be 
written in terms of its real and imaginary parts: 

W fg (a,T) = CoW fg {a,T) - iQuadW fg {a,T). (3) 

It can be shown that the following equation holds if f(t),g(i) £ £ 2 (3fJ) 

/ + OC Z' + OC /■ + OO 

f(t)g(t)dt=l/c* / CoW fg {a,T)dTda 7 (4) 
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where is a constant depending on the choice of the wavelet mother. 



3 Cross wavelet coherence functions 

The highly redundant information from a multiscale wavelet analysis of time 
series must be reduced by means of suitable selective procedures and quantities, 
in order to extract the main features correlated to an essentially intermittent 
dynamics. In this study, we analysed and compared the properties of two com- 
plementary wavelet local correlation coefficcnts which are able to well evidence 
peculiar and anomalous local events associated to the vortex dynamics. More 
precisely, given two signals f(t) and g(t), we refer to the so-called Wavelet Local 
Correlation Coefficent (Burcsti et. al |lj]), defined as: 

WLCC { a,r)= . w f 0W ^ .. . (5) 
I W f (a,T) || Wg(a,T) | 

This quantity is essentially a measure of the phase coherence of the signals. Here 
we introduce the Cross Wavelet Coherence Function (CWCF) defined as: 

I Wf(a,T) | 4 + | Wg{a,T) | 4 

which is essentially a measure of the intensity coherence of the signals. Using 
the polar coordinates we can write the wavelet transforms of Wf(a, t), W g (a, r) 
and Wf g (a, r) as: 

W f (a,r)=p f e^ W g {a, T ) = p g e«» (7) 

W fg (a,r) = p fPg e< e '-°t\ (8) 
and the Cross Wavelet Coherence Function can be written also as: 

CWCF(a,r) = -f^j. (9) 

Pf "I" Pg 

It is easy to observe the two basic properties of the function (6): 

CWCF{a, T ) = 0^p f = or p g = (10) 
0<CWCF <1 V(»,t, (11) 
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4 Numerical simulation 

Wc considered a laboratory-scale axisymmctric flame of methane-air in a non 
confined bluff-body configuration. More precisely, the burner consists of a 5.4 
mm diameter methane jet located in the center of a 50 mm diameter cylinder. 
Air is supplied through a f 00 mm outer diameter coaxial jet around the 50 mm 
diameter bluff-body. The Reynolds number of the central jet is 7000 (methane 
velocity =21 m/s) whereas the Reynolds number of the coaxial jet is 80000 (air 
velocity =25 m/s). This is a challenging test case for all the turbulence models, 
as well documented in the ERCOFTAC report (Chatou, 1994) ||. Moreover, due 
to the highly intermittent, unsteady dynamics involved and the high turbulence 
level, especially for the reactive case, the Large Eddy Simulation (LES) appears 
as the most adequate numerical approach (Sello et. al Q). 

5 Results and discussion 

In this analysis we are mainly interested to relations existing between evolution 
of turbulence and mixing, for the reactive case. Previous DNS simulations on 
coaxial jets at different Reynolds numbers, show the ability of the wavelet cross- 
correlation analysis to better investigate the relations between mixing process 
and the dynamics of vorticity (Salvetti et. al Q|). Thus, the signals analysed here 
arc velocity fluctuations (for Reynolds stress contributions) and mixture fraction 
fluctuations (for mixing evolution) from LES. As an example, Figure 1 shows the 
wavelet co-spectrum maps for a significant time interval in the pseudo-stationary 
regime of motion. The main contributions to the Reynolds stress are evidenced 
by high intensity correlations (red) and anti-correlations (blue) regions, which 
evolve intermittently. The dominant frequencies involved are located around 
130 Hz. For the mechanisms responsable of the evolution of mixing, we note 
that the same regions of high Reynolds stress correspond to high correlation, 
or cooperation, between velocity and mixture fraction fluctuations, suggesting 
that, at the selected location, the same events of stretching and tilting of the 
vorticity layer, drive both Reynolds stress and mixing evolutions. Note that the 
large high value region located at low frequencies in the right map is statistically 
not significant if we assume a proper red noise background spectrum. To better 
investigate the role of the high correlation regions, we performed a cross section 
in the wavelet map at the frequency 160 Hz. Figure 2 (left) shows the time 
behaviour of the coherence functions WLCC, eq.(5), and CWCF, eq.(6). Here 
the phase and intensity coherence of signals are almost equivalent, but we can 
clearly point out an important anomalous event occurred at around t=0.19 s, 
corresponding to a loss of both intensity and phase coherence, followed by a 
change of the correlation sign. The link between this event and the dynamics 
of vorticity is evidenced by Figure 2 (right), which displays the wavelet map of 
the related vorticity signal. The higher frequency significant regions (~ 730 Hz) 
result strongly intermittent, with a bifurcation to lower and higher values than 
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average, followed by a drop of activity, in phase with the anomalous event. 



Figure 1: Cross- Wavelet co-spectrum maps for axial and radial velocity fluctu- 
ations (left) and for axial velocity and mixture fraction fluctuations (right) at a 
given spatial point near the edge of the central jet. 



Figure 2: Coherence functions for axial velocity and mixture fraction fluctuations 
(left) and wavelet map of vorticity time series (right). 

These few examples support the usefulness of the cross- wavelet analysis ap- 
proach to better investigate turbulent mixing processes in real systems. 
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